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Abstract 

We study the conductance through finite Aharonov-Bohm rings of interacting electrons weakly coupled to non- 
interacting leads at two arbitrary sites. This model can describe an array of quantum dots with a large charging 
energy compared to the interdot overlap. As a consequence of the spin-charge separation, which occurs in these highly 
correlated systems, the transmittance is shown to present pronounced dips for particular values of the magnetic flux 
piercing the ring. We analyze this effect by numerical and analytical means and show that the zero-temperature 
equilibrium conductance in fact presents these striking features which could be observed experimentally. 
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1. Introduction 

One of the challenges of nanotechnology is the pos- 
sibility of fabricating new artificial structures with 
tailored properties. For example, the Kondo effect 
was achieved in a system consisting of one quantum 
dot connected to leads[f ,2,3]; systems of a few QD's 
have been proposed theoretically as realizations of the 
two-channel Kondo model [4,5], the ionic Hubbard 
model, [6] and the double exchange mechanism. [7] 
Also, the correlation-driven metal-insulator transition 
has been studied in a chain of quantum dots. [8] 

Another interesting phenomenon in strongly corre- 
lated systems is what is known as charge-spin separa- 
tion. It is well known that strong correlations in one 
dimension invalidate the Fermi liquid conventional de- 
scription of electrons. In particular, correlations can 
lead to the fractionalization of the electron into charge 
and spin degrees of freedom. [9] This separation is an 
asymptotic low-energy property in an infinite chain. 
However, exact Bethe ansatz results for the Hubbard 
model in the limit of infinite Coulomb repulsion U 
show that the wave function factorizes into a charge 
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and a spin part for any size of the system. [10] There 
have been several experiments reporting indirect in- 
dications of charge-spin separation [11,12,13], and it 
could also be potentially observed in systems such as 
cuprate chains, ladder compounds, [14] and carbon 
nanotubes. [15] 

Several theoretical approaches tackled this phe- 
nomenon in the ring geometry. For example, the real- 
time evolution of electronic wave packets in Hubbard 
rings has shown a splitting in the dispersion of the 
spin and charge densities as a consequence of the dif- 
ferent charge and spin velocities. [16,17] Pseudospin- 
charge separation has also been studied in quasi-one- 
dimensional quantum gases of fermionic atoms. [18,19] 
Other theoretical approaches concern the transmit- 
tance through Aharonov-Bohm rings modeled by 
Tomonaga-Luttinger liquid or other correlated Hamil- 
tonians like the Hubbard or t — J models. [20,21,22,23] 
In these systems, noticeable dips are observed in the 
conductance for fractional values of the flux which can 
be attributed to destructive interference between de- 
generate states as we explain below. We have recently 
discussed the extension of these results to ladders of 
two legs as a first step to higher dimensions. [24] 

In this work, we analyze the origin of the dips in 
the transmittance as a function of applied flux in finite 
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Fig. 1. Schematic representation of the systems studied numer- 
ically, (a) L = 6, (b) L — 7. In both cases M — 3. 

rings described by the t — J model. We discuss the 
conditions for which the intensity of the lowest-lying 
peak in the zero-temperature equilibrium conductance 
as a function of the gate voltage presents characteristic 
dips for certain flux values. These will be shown to be 
a consequence of spin-charge separation. 



2. Model Hamiltonian 

The basic model is depicted in Fig. 1. We have con- 
sidered a ring of L sites, weakly connected to non- 
interacting leads at sites and M. 

The Hamiltonian reads: 



H — -Wring + Pleads + blinks- 



(1) 



The first term describes the isolated ring, with on-site 
energy given by a gate voltage V g , and hoppings mod- 
ified by the phase exp(i<j)/L) due to the circulation of 
the vector potential. For most of the results of this pa- 
per we use the t — J model to describe the ring, 

#ring = -eVg c\a C i° ~ 1 ( C I+1 CT C i° + H ' C -) 

icy 

+jWs i .s i+ i-iy (2) 

where (j> = 27r<I>/ < l? u , /$o = h/2e is the flux quantum, 
Si is the spin operator at site i and double occupancy 
is not allowed at any site of the ring. The second term 
corresponds to two tight-binding semi-infinite chains 
for the left and right leads 



3. Conductance 

When the ground state of the isolated ring is non- 
degenerate, and the coupling t' between the leads and 
ring is weak, the equilibrium conductance at zero tem- 
perature can be expressed to second order in t' in terms 
of the retarded Green's function for the isolated ring 
between sites i and j: Gfj(co). [6,20] For an incident 
particle with energy lo = — 2t cos k and momentum ±fe, 
the transmittance reads 



T(<j,V g ,4>) = 



At 2 sin 2 ifc|t(w)| 2 



[w - e(w) + te ik ] 2 - \i 2 (u) 



(5) 



where e(u>) = t' 2 GoQ{oj) and t(u>) = t' 2 G^ M (oj) repre- 
sent a correction to the on-site energy at the extremes 
of the leads and an effective hopping between them re- 
spectively. 

This equation is in fact exact for a non-interacting 
system, however, it loses validity for an odd number of 
electrons, where the ground state is Kramers degener- 
ate. In this case the method misses completely the in- 
teresting physics arising from the Kondo effect. [6,25] 
However, the Kondo temperature is small compared 
to the other energy scales of the system and we can 
assume safely that the Kondo effect is destroyed by a 
small temperature. This approach is justified for small 
enough t' since the characteristic energy of this Kondo 
effect decreases exponentially. 

The conductance is G = (ne 2 /h)T(n,V g ,4>), where 
n — 1 or 2 depending if the spin degeneracy is broken 
or not, [6] and p is the Fermi level, which we set as zero 
(half-filled leads). When the gate voltage V g is varied 
a peak in the conductance is obtained when there is 
a degeneracy in the ground state of the ring for two 
consecutive number of particles: E g (N + 1) = E g (N), 
where E g (N) is the ground state energy of H T i ng with 
iV electrons. Without loss of generality, we assume that 
we start with N + 1 electrons in the ring and apply a 
negative gate voltage in such a way that a peak in the 
conductance is obtained at a critical value V g when the 
number of electrons in the ring changes from iV + 1 to 
iV electrons. 



4. Analytical and numerical results 
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The third term in Eq. (1) describes the coupling of the 
left (right) lead with site (M) of the ring 



-Hiinks = -t' y^oli.gCpg + a\ a c M o + H.C.). 



(4) 



For J = the model is equivalent to the Hubbard 
model with infinite on-site repulsion U, for which the 
wave function can be factorized into a spin and a charge 
part, evidencing charge-spin separation. [10,22,26] For 
each spin state, the system can be mapped into a spin- 
less model with an effective flux which depends on the 
total spin. For a system of iV particles one can construct 
spin- wave functions with wave vectors k s — 2irns/N, 



2 



where the integer n s characterizes the spin wave func- 
tion. The total energy and momentum (in an appro- 
priate gauge) of any state of the ring have simple ex- 
pressions: 



E= -2t^cos(fci 



2nni + (j> e s 



K = Y^ k i = I 27r K + n s ) + N(t>] /L, 



i , ± 2tt 
off = <P + k s = (p + —n B , 



(6) 

(7) 
(8) 



where the integers ni characterize the charge part of 
the wave function and n c = Y)ni. 

The values of the flux for which dips or reduced 
conductances are expected, correspond to some partic- 
ular crossings of the energy levels of N electrons. One 
can see this from the general form of the Green's func- 
tions Goj(uj) entering the transmittance [Eq. (5)] when 
a particle is destroyed. Using the Lehman's represen- 
tation, the relevant part of the Green's function is: 



UJ + E e - Eg 



0) 



Noting that the ring has translational symmetry and 
naming K v the wave vector of the state \v), one obtains: 



0jy 1 ^ tu + E e -E g 



(10) 



At certain flux values, two states of N electrons, 
\e) and |e'), become degenerate. Assuming that the 
corresponding matrix elements entering Eq. (10) are 
nonzero, only these two states, contribute significantly 
to the Green's function Goj(aj) at the Fermi energy 
{u) = fj. = 0) when V g (which displaces all E e rigidly 
with respect to E g ) is tuned in such a way that E e = 
E e , ~ E g . Denoting by /3 = exp[iM(K e > - K e )] the 
relative phase between the two intervening degenerate 
states we see that if f3 7^ 1, the transmittance, which is 
proportional to \G^ M (u})\ 2 , (Eq. (5)), is reduced near 
the crossing (note that this results is independent of 
any specific model). 

We can predict the positions of the dips in the trans- 
mittance without having to resort to the calculation 
of the matrix elements entering the Greens functions 
Eq. (10). The crossings of energy levels at low energies 
take place at cj> = —n(n s + n' s )/N. [23] When n s + n' g 
is odd (even) the relative phase /3 = exp[iL(K e * — 
K e )/2] = exp [i(n' s — n s )] = —1 (1) and there is (there 
is not) a dip in the integrated transmittance. There- 
fore, the positions of the dips are located at 
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Fig. 2. Intensity of the first peak in the transmittance If p , as 
a function of applied magnetic flux, for a ring of 6 sites and 
(a) iV + 1 — 5 and (b) N + 1 — 6 electrons in the ground state, 
three different configurations M = 1, 2 and 3, t' — Q.3t and 
.7 = 0.001*. 



with n integer. These are also the positions where cross- 
ings in the (experimentally accessible) ground state for 
N particles take place (n' s — n s = ±1). 

To check these predictions we have performed nu- 
merical calculations for the transmittance, obtained by 
diagonalizing the ring using Davidson's method [27]. 
Once the Green's functions were obtained, they were 
replaced in Eq. (5), to obtain the transmittance. The 
systems studied are represented in Fig. 1. In contrast 
to previous work, [20,21,22,24] we concentrate on the 
first peak in the transmittance as the gate voltage is 
decreased since this is the feature which is experimen- 
tally accessible at equilibrium and low temperatures. 

In Fig. 2 we show numerical results for a ring with 
L — 6 sites and N = 4 and 5 particles in the inter- 
mediate state for the three non-equivalent configura- 
tions corresponding to M = 1,2 and 3. To quantify the 
relative intensity of the conductance we integrate the 
transmittance given by Eq. (5) over a window of gate 
voltage V g of width 0.002t centered around the degen- 
eracy point between the ground state for N + 1 and 
iV electrons. This corresponds to the intensity of the 
first observable peak in the transmittance as the gate 
voltage is lowered. As the curve is symmetric under 
change of sign of <f> we show only the interval < cj) < 
7T. For ,7^0 and N = 4 the dips appear as expected 
(Eq. (11)) at </> d /TT = 0.25 and 0.75. For N = 5 the 
dips should occur at 4>d/^ = 0.2, 0.6 and 1. However, 
near 0.6 the total spin in the ground state for 5 elec- 
trons changes from S = 1/2 to S — 3/2, which is not 
accessible by destroying an electron in the 6-electron 
singlet ground state. Therefore, the transmittance van- 
ishes for 0.6 < < 7r. As a consequence, in the inter- 
val shown there is only one dip present. In this figure 
one can see the effect of the different source-drain con- 
figurations: the lowest conductance is achieved for the 
symmetric case (M — 3 for L — 6) where j3 = — 1 and 
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Fig. 3. Intensity of the first peak in the transmittancc If p , as 
a function of applied magnetic flux, for a ring of 7 sites and 
several fillings. Here M = 3, t' = 0.3t and J = O.OOlt. 



both degenerate levels interfere destructively. For the 
other cases the interference is less destructive and a 
less pronounced dip is obtained. 

It is also interesting to analyze the conductance for 
different particle numbers to study the evolution of the 
dips as predicted in Eq. (11). This is shown in Fig. 3. 
The number of particles in the intermediate situation, 
N, is shown for each case. Here we see that the minima 
do in fact occur for the fluxes predicted by that relation 
very accurately. The abrupt step obtained for N = 2 
corresponds, again, to a forbidden transition to a large 
total spin state from a singlet ground state as explained 
above. 



5. Conclusions 

We have presented results for the zero-temperature 
equilibrium conductance through finite rings described 
by the t — J model threaded by a magnetic flux, weakly 
coupled to conducting leads. At particular values of the 
flux we find dips or reductions of the transmittance, 
which are due to negative interferences between de- 
generate levels. This can be understood by analyzing 
the extremely interacting case for J = 0, where exact 
results are available. The position of the dips reflect 
the particular features of the spectrum in this limit, in 
which the charge and spin degrees of freedom are sep- 
arated at all energies. For finite J the positions of the 
dips change and some additional dips can also appear 
in a manner that is difficult to predict and which is not 
yet fully understood. 

The negative interference depends on the source- 
drain configuration. It is more marked if the leads are 
connected at angles near 180 degrees. These results are 
confirmed by our numerical calculations. 
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